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Abstract — In this paper, a novel technique is applied to detect 
fault in the transmission line using wavelet transform. Three 
phase currents are monitored at both ends of the transmission 
line using global positioning system synchronizing clock. 
Wavelet transform, which is very fast and sensitive to noise, is 
used to extract transients in the line currents for fault 
detection. Fault index is calculated based on the sum of local 
and remote end detail coefficients and compared with 
threshold value to detect the fault. Proposed technique is 
tested for various faults and fault inception angles. Simulation 
results are presented showing the selection of proper 
threshold value for fault detection. 

Index Terms — Wavelet transform, transmission line faults, 
power system protection, fault transients, multiresolution 
analysis 

I. Introduction 

Transmission lines exposed to various faults are designed 
with distance protection scheme. The main principle of this 
technique is to calculate the impedance between the relay 
and fault points; the apparent impedance is then compared 
with the relay trip characteristic to ascertain whether it is an 
internal or external fault. Microprocessor based distance 
relays have been broadly applied to transmission lines 
protection. The calculation of line impedance has been carried 
out in difference equations [1], [2] or using phasor obtained 
by Fourier algorithms [3], [4]. To set a distance relay, it is 
necessary to determine the reach of its setting. A distance 
relay does not aim to locate a fault but instead its objective is 
to determine whether a fault is within or outside the intended 
protection zone. Methods exploiting travelling wave 
components to achieve ultra high speed protection have also 
been proposed [5], [6], while their performance needs further 
evaluation and, where necessary, improvement [7] . Wavelet 
analysis as a powerful tool for signal processing can be 
applied to effectively overcome difficulties of the travelling 
wave protection techniques [8]-[l 1]. By virtue of the variable 
length data window of the wavelet transform (WT), it possible 
to detect and characterize closed singularities of 
nonstationary signals. The wavelet transform can also be 
employed for preprocessing of input data to improve the 
performance of artificial-neural-network (ANN)-based 
algorithms [11]. In [12] principal component analysis 
technique is proposed to identify the dominant pattern of 
the signal preprocessed by the wavelet transform. This way, 
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the extracted information of the voltage and current signals 
by wavelet transform at different scales is utilized. This 
technique improves the capability of traveling-wave 
detection especially for the faint and close-in faults. The 
proposed position protection algorithm is based on polarity, 
magnitude, and time interval between the arriving waves, 
and can rapidly identify most of the faults occurring inside 
the protected zone. In this paper, wavelet based 
multiresolution analysis is used for detection of fault based 
on detail- 1 coefficients of current signals at both the ends of 
line. Transients are introduced due to fault in the currents of 
faulted phase. These high frequency components are easily 
extracted by wavelet analysis. A fault index is then calculated 
by combining the detail- 1 coefficients of three phase currents 
monitored at both ends of the transmission line using global 
positioning system synchronizing clock. The details of the 
proposed scheme are described and simulation results are 
presented for detection of various faults at different fault 
distance and inception angles. Results show that a common 
threshold fault index is used to detect the fault in the line. 

II. Wavelet Transforms 

Wavelet transform is useful in detecting and extracting 
disturbance features of various types of electric power quality 
disturbances because it is sensitive to signal irregularities. 
Wavelet analysis expands functions not in terms of 
trigonometric polynomials but in terms of wavelets, which 
are generated in the form of translations and dilations of a 
fixed function called the mother wavelet. A mother wavelet is 
a function that oscillates, has finite energy and zero mean 
value. Compared with Fourier transform, wavelet can obtain 
both time and frequency information of signal, while only 
frequency information can be obtained by Fourier transform. 
The signal can be represented in terms of both the scaling 
and wavelet function [8] as follows: 
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Cj is the J level scaling coefficient 
d is the j level wavelet coefficient 
ij>(t) is scaling function 
^(t) is wavelet function 
J is the highest level of wavelet transform 
t is time 
Each wavelet is created by scaling and translation operations 
of mother wavelet. Wavelet theory is expressed by 



-fcACEEE 



ACEEE Int. J. on Electrical and Power Engineering, Vol. 02, No. 02, August 201 1 



continuous wavelet transformation as 

00 * 

CWT ¥ x(a,b) = W x (a,b)= \ x(t)y ab (t)dt (2 ) 

— 00 

where W „ h (0 =1 a I 1 ' 2 y/ (^— ) , 

a (scale) and b (translation) are real numbers. 
For discrete-time systems, the discretization process leads 
to the time discrete wavelet series as 



DWT x(m,n)= f x(t)y/ (t)dt 
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In Multi resolution analysis (MRA), wavelet functions and 
scaling functions are used as building blocks to decompose 
and construct the signal at different resolution levels. The 
goal of MRA is to develop representations of a signal at 
various levels of resolution. MRA is composed of 2 filters in 
each level which are low pass filter and high pass filter. MRA 
can be shown as in Fig. 1. Wavelet coefficients obtained 
from high pass filter are called detail coefficients and 
coefficients at low pass filter are called approximation 
coefficients. Transients generated during faults are of high 
frequency components which can easily be extracted from 
level 1 wavelet coefficients. 

III. Fault Protection Algorithm 

In this paper, wavelet analysis is used to detect the fault 
in the transmission line. Fig. 2 shows system configuration 
used for analysis. A 100 kms transmission line of rating 138kV, 
100MW, 60Hz is fed from both ends. Synchronized sampling 
of three phase currents and voltages is carried out at both 
the ends with the help of a GPS satellite. Wavelet analysis of 
these sampled signals is performed and level- 1 detail 
coefficients are used to calculate fault index and detect fault. 
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Figure 1. Four level multi resolution signal decomposition of a 
signal. 



a 



rinriadukii l.ln- 



I V*l 



JO 



WTbaied 

fault detector 



r 



WT based 

f»uit dttecior 



Figure 2. Transmission line system with FACTS device 
implemented with proposed scheme 
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Figure 3. Current waveforms during ABG fault at 20% line length 

from bus 1 (a). Three-phase currents at Bus-1 and (b). Three-phase 

currents at Bus-2 

Two-phase to ground fault is considered here for fault 
index calculation. Figs. 3(a) and (b) show the three-phase 
currents at bus-1 and bus-2, respectively for ABG fault. 
Currents in phase A and B increase during fault. These 
sampled currents are then analyzed with Debauches' wavelet 
to obtain level- 1 detailed coefficients (CD1) over a moving 
window of half-cycle length. Level- 1 detail coefficients for 
phase-A, B and C currents (CDl_ial, CDl_ibl, CDl_icl) at 
bus-1 are shown in Figs. 4(a), (b) and (c) respectively. Level- 
1 detail coefficients for phase-A, B and C currents (CDl_ia2, 
CDl_ib2, CDl_ic2) at bus-2 are shown in Figs. 4(d), (e) and 
(f) respectively. These coefficients indicate high frequency 
components in the currents due to fault. Currents in faulted 
phase indicate a large magnitude of coefficients which are 
used to detect fault. 
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Figure 4. Wavelet analysis of three-phase currents, (a), (b), (c). Level-1 detailed coefficients of bus-1 phase-A, phase-B current and phase-C 

currents respectively, (d), (e), (f)- Level-1 detailed coefficients of bus-2 phase-A, phase-B and phase-C currents respectively, (g), (h), (i). 

Effective detail coefficients of phase-A, phase-B current and phase-C currents respectively 



The detail coefficients received from the remote bus (bus-2) 
are added to the local bus (bus- 1 ) detail coefficients. Variation 
of these effective detail coefficients (CDl_ia, CDl_ib, CDl_ic) 
are shown in Figs. 4(g)-(i). 
Fault index is then calculated in each phase as 

Ifl=T\DI E \ (4) 

Variation of fault index in each of the three-phases for ABG 
fault is shown in Fig. 5. It is clear from the figure that fault 
index is large in faulted phases as compared to that in healthy 
phase. A threshold value (I m ) is then set such that I m < I of 
faulted phase and I m > I of healthy phase. This algorithm is 
tested for various faults and a common threshold value is 
selected. 

IV Simulation Results 

The proposed algorithm is tested for all types of faults, 
considering variations in fault locations and fault inception 
angles in the range of 0-1 80° for the simulation circuit of the 
system shown in Fig. 6. Simulation is performed for an LLG 
fault at every 20kms of the line and at various fault inception 
angles. Fig. 7 shows variation of fault index, / for ABG fault 
applied at every 20kms distance on the line and fault inception 
angles of 0°, 30°, 60°, 90°, 150° and 180°. It is clearly observed 
from the figure that fault index is large in faulted phases A 
and B compared to fault index in healthy phase C. This shows 
that extracting transients from the fault currents using 
wavelets helps to detect fault and faulted phase. To compare 
the fault index for various types of faults, simulation is 
performed for all the faults. Fig. 8 shows the variation of fault 
index for AG BG ABG AB, ABCG and ABC faults at different 
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distances and for fault inception angle at 30°. It can be seen 
that transients in faulted phases are high compared to healthy 
phase. From Figs. 7 and 8 it is possible to set a threshold 
value between the fault index for fault and healthy phases. 
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Figure 5. Variation of Fault Index in due to ABG fault 
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Figure 6. Simulation circuit of the transmission system used for 
testing the proposed algorithm 
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x-axis: Fault distance (kms) 
y-axis: Fault index 



* Phase A fault index, Phase B fault index, + Phase C fault index 



ABG 180deg 
O * 



100 

(f) 



200 



Figure 7. Variation of fault index for ABG fault 
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Variation of fault index connected and fault at 140kms 
distance from busl for various faults 



Fault index is calculated for case studies of all possible faults 
at various points on the line and at various inception angles. 
These values are tabulated in Table I The values listed are 
minimum fault index possible for faulty phase and maximum 
fault index possible for healthy phase. For example 
considering ABG fault from Table 1, minimum fault index in 
phases A and B are 2.5397 and 2.0021 respectively and 
maximum fault index in phase C is 0.4557. Fault index calculated 
from transients show that index values are larger for faulty 
phase currents compared to index of healthy phase current. 
The proposed algorithm is therefore having an edge over 
fundamental impedance calculation method which is normally 
employed with distance relay. From Table 1 minimum fault 
index in phase A (/,, ) is 1.3589 which is obtained 

■*■ s jia_min' 

considering for all possible faults with phase A involved i.e., 
AG fault, ABG fault, CAG fault, AB fault, CA fault, ABCG 
fault and ABC fault. Maximum fault index in phase A (/„ ) 
is 0.5794 which is obtained for faults with healthy phase A 



i.e., BG fault, CG fault, BCG fault and BC fault. Similarly 
minimum and maximum fault index for Phase B and C (/ 
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I ) are identified. Proper threshold index 
(I TH ) is then set between minimum fault index (I fl .) and 
maximum fault index (I ) of all faults as 



flju 
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where /„ . = min(/„ . , /„, . , /„ . ) and 

jljmn v flujnuv flbjnm flcjmn' 

/,, =max(/„ ,/„, ,/„ ) 

jljnux v flajnux flbjnux flc_max / 

From this case study, /„ . = 0.7393 and /_ = 0.5794 (from 
Table 1). A common threshold index / = 0.65 may be set 
between /„ and /„ . Fault in a phase is detected when 

fl_imn jl_max x 

the fault index exceeds 0.65. The feature of this algorithm is 
that transients are easily extractable using wavelet analysis 
and provides a common setting of threshold index. 

Table I. Fault Indices for Various Faults 



Fault 


Phass-A Index 


Phass-B Index 


Phase-C Index 


Mm 


Max. 


Mm 


Max. 


Min. 


Max. 


AG 


IJ589* 


- 


- 


0.5715" 


- 


0.555-e- 


BG 


- 


0.49SS 


0.7690* 


- 


- 


0.4337 


CG 


- 


0.5794" 


- 


0.5577 


0.7393* 


- 


ABG 


2.5397 


- 


:.cc:: 


- 


- 


0.4557 


BCG 


- 


0.5179 


11915 


- 


0S0S4 


- 


CAG 


LS323 


- 


- 


0.4447 


16221 


- 


AB 


2.335C 


- 


2.3532 


- 


- 


0.0253 


BC 


- 


0.0543 


3212 


- 


03315 


- 


CA 


1.S047 


- 


- 


0.0405 


13366 


- 


ABCG 


2.7743 


- 


1.7337 


- 


1.2769 


- 


ABC 


2.7743 


■ 


1.7337 


- 


1.2769 


■ 



: indicates the lowest fault index when the phase is faulted 
indicates the largest fault index when the phase is healthy 
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V. Conclusions 

An algorithm for detection of faults in the transmission 
line is proposed based on transient component extraction of 
currents. Fault index calculated from detail- 1 coefficients of 
currents are clearly separable for fault and healthy phases. 
This provides an easy selection of threshold to detect fault 
in the line. Simulation results presented for various faults at 
different distances and inception angles show that a common 
threshold is set to detect a fault in the transmission line. 
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